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Abstract 

The B(E2:0 1 ~ — ► 2^) values of the Ba isotopes (Z = 56) exhibit a sharp 
increase in deformation as the neutron numbers approach the mid-shell value 
of N = 66. This behavior is anomalous because the 2f level energies are 
very similar to those of the neighboring isotopes. By means of the axially- 
symmetric deformed Woods-Saxon ( WS) hamiltonian plus the BCS method, 
we investigated the systematics of B(E2) of the Ba isotopes. We showed that 
15% of the B(E2) values at N = 66 was due to the level crossing, occurring 
at the deformation /?ws ~ 0-3, between the proton orbits originating from 
the orbits 9? = l/2~(Ml/2) and 9/2+(st9/2) at zero deformation. The latter 
of these two was an intruder orbit originating from below the energy gap at 
Z = 50, rising higher in energy with the deformation and intruding the Z = 
50 — 82 shell. These two orbits have the largest magnitude of the quadrupole 
moment with a different sign among the orbits near and below the Fermi 
surface. Occupancy and non-occupancy of these orbits by protons thus affect 
B(E2:0^ -> 2f) significantly. 

PACS Numbers: 21.10.-k; 21.10.Ky; 21.60.Ev; 21.60.Fw; 23.20.-g 
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I. INTRODUCTION 

The deformation of the nuclear ground state is a fundamental quantity that can be 
deduced from the B(E2) value for the first 2 + state (J 71 = 2f). From the systematics of 
the B(E2:0 ] f — > 2^ ) values, we could infer how the deformation of the ground state changes 
from nucleus to nucleus. Over the past decade, the B(E2) values of the Xe-Ba-Ce nuclei 
have been measured extensively. See for Xe, @ for Ba, || for Ce. The accumulated data 
has enabled us to compare the data with several theoretical results to test these theories [|J . 

Focusing ourselves on the Ba (Z=56) isotopes, we will discuss the systematics of the 
B(E2:0^ -> 2f) values. 



II. GRODZINS RULE 

The B(E2) values of Ba []5J as well as the Xe (Z=54) isotopes || increase sharply as the 
neutron number N approaches N = 66 |4]|| as can be seen in Fig. [I]. The solid squares 
in Fig. m present the experimental data of B(E2:0^ — > 2+), which are compared in Fig. [I] 
with the theoretical results estimated using the Grodzins rule [§] and obtained using the 
proton-neutron interacting boson model (IBM-2) [|7|, as we will explain in detail below. 

The Grodzins rule expresses an empirical relation between Ex(2+) and B(E2): 

B(E2 : 0+ -> 2 + ) = 1.63 x 10~ 2 Z 2 /(£ ■ A){ebf (1) 

with E being Ex(2^) (MeV), Z the proton number, and A the mass number. In deriving 
Eq. (0), we used Eq. (4) of Ref. ||. The open triangles in Fig. Q present B(E2:0^ — ► 2*) 
derived from Eq. ([!]). The results obtained by this rule are in surprisingly good agreement 
with the data for the less deformed Ba isotopes with iV > 70 but underestimate B(E2) by 
20% at iV = 66. The similar empirical relation in [|J also underestimates B(E2:0 ] f — > 2f) at 
N = 66. The empirical rules suggest that the sharp increase of B(E2) is anomalous because 
Ex(2i ) at iV = 66 is almost the same as those of the neighboring isotopes. We will next 
examine this anomaly in B(E2:0^ — > 2f ) within the framework of IBM ||. 



III. IBM-2 

The interacting boson model (IBM) || has been widely used for describing the 
quadrupole collective states of the medium heavy nuclei. The microscopic foundation of 
IBM has been given within the proton-neutron IBM (IBM-2) |§. The building blocks of 
IBM-2 are s T and d T bosons (r = it, v) which are considered to be approximations to the pro- 
ton (neutron) pairs with spin-parity + and 2 + . The boson images of the fermion operators 



are given in terms of the OAI mapping [10 



Using IBM-2, Otsuka, Pan and Arima |7j investigated the Xe-Ba-Ce isotopes concentrat- 
ing on the systematics of B(E2:0 ] f — > 2f) as a function of the neutron number. According to 



the procedures given in flTfl, we did the IBM-2 calculations using the code NPBOS [11 1 with 
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the parameters taken from Ref. pj which were derived from the microscopic calculations 



through the method of OAI mapping [10]. The boson image of the fermion E2 operator is 
given by pj 



(2) 



where e° (a = tt, u) are constants and 



(2) 



(3) 



with 



\ (0,-1) 



(4) 



X<x = q _ 2 " K * (5) 

where is the degeneracy of the major shell and = 16 for the Ba isotopes. We used the 
following parameters taken from J7J: Xa = ±0.672 (— for the particles and + for the holes), 
e° = 0.154 (eb), e° = 0.110 (eb). Using this E2 operator, we calculated B(E2:0| -> 2f) for 
the Ba isotopes, and compared them with the experimental data |2| in Fig. [I]. The solid line 
represents the theoretical results which were identical to the results in |7|] . 

Quantities K a and Xa/Xa represent the effects of the Pauli principle, and are referred to 
as the Pauli factors in 0. These Pauli factors become unity when there is no Pauli effect, 
i.e., Q a — > oo 0. The dashed line in Fig. |l] is the B(E2) values calculated without the Pauli 
factor n a in §, i.e., with K a = 1 and Xo = xlPa ~ 2N a )/{Vt a - 2). 

In Fig. |l], we used the latest experimental B(E2) data, some of which || were measured 
after the publication of [0. We should mention again that, besides the experimental data, 
the IBM-2 results in Fig. g are identical to Fig. 4 of because we have taken the same 
parameters as was used in 0. 

It is clearly seen from Fig. [I] that the IBM-2 predictions with the Pauli factor are in 
good agreement with the data for N > 76. However, for 70 < N < 74, IBM-2 without the 
Pauli factor gives B(E2) values closer to the data than IBM-2 with the Pauli factor. For 
iV = 66,68, the B(E2) values of IBM-2 with the Pauli factor saturate in the mid-shell in 
disagreement with the experimental trends. At iV = 66, the IBM-2 calculations without the 
Pauli factor yield nearly the same magnitude of B(E2:0 ] f — > 2f) as the Grodzins rule. Thus, 
at iV = 66, not only the Grodzins rule but also IBM-2 underestimate the B(E2) values. 

The Xe isotopes also showed a sharp increase in the B(E2) values at the mid-shell [|TJ. 
Raman et al. pointed out in Ref. M that the single-shell model such as the single-shell 
Nilsson model failed in reproducing the B(E2) values of the Xe isotopes at iV ~ 64-66, 
whereas several multi-shell models correctly predicted these values. They showed in || that 
the rapid rise of the Q = 9/2 (g9/2) oblate orbit for protons as a function of deformation and 
its intrusion into the 50-82 shell suggests a mechanism for increasing the B(E2:0 ] f — > 2^) 
value as a result of partial emptying of this orbit due to pairing. However, the detailed 



3 



manner by which the mid-shell Xenon isotopes acquire significant deformation is still unclear, 
as pointed out in |4j]. 

We aim in this paper at clarifying the origin of the sharp increase of B(E2:0^ — > 2f) 
observed in the mid-shell Ba isotopes. To achieve this objective, we used the axially sym- 
metric deformed Woods-Saxon (WS) hamiltonian plus the BCS method. All orbits with 
eigenenergies lower than the barrier height [12| were adopted in the calculations. Therefore, 
we can get rid of the drawbacks || of the single-shell models such as the single-shell Nilsson 
model. In contrast to the several theoretical tools examined in ||, we do not minimize 
a total hamiltonian to fix the intrinsic state of the ground band. Instead, as a function of 
deformation parameter /? W s of the WS potential, we calculated intrinsic quadrupole moment 
Qq. Then, we determined /?ws directly from the experimental Qq values, by which we could 
discuss the variation of the intrinsic configuration of the ground band with the deformation. 
We will show that 15% of B(E2:0+ — > 2f) for the Ba isotope with N = 66 is due to the level 
crossing between the Q = 9/2(^9/2) and Q = l/2(/ill/2) orbits, but the effect of the level 
crossing on B(E2) fades out rapidly as N goes away from the mid-shell value of N = 66. 
This suggests that the sharp increase in B(E2) at the mid-shell is due to the level crossing 
between these orbits. 



IV. WOODS-SAXON HAMILTONIAN PLUS BCS CALCULATIONS 

We assume in this paper that individual nucleons move in the axially-symmetric 
quadrupole deformed Woods-Saxon potential and interact with each other by the monopole 
pairing force. We take the same form of the Woods-Saxon (WS) potential as was used in 
fl~2|| . We assume Y 2< o deformation only, and denote the deformation parameter of the WS 
potential by /3ws as was already mentioned. 

We should mention here the 7-softness or triaxiality of the heavier Ba isotopes with 
iV < 82. The Ba isotopes exhibit the 7-soft or 0(6)-like level scheme for the heavier isotopes 
with iV < 82 |13|,14], and change gradually from the 7-soft to axially symmetric nuclei as 



N approaches the mid-shell value of iV = 66 |L4]]. Through this paper, we assumed the 
axially symmetric WS potential for simplicity. This assumption may not be so valid for the 
heavier Ba isotopes because of the 7-softness. However, our major interest lies in the rapid 
increase in B(E2) of the Ba isotopes with iV ~ 66. For these nuclei, the axial symmetry is 
a reasonable approximation. 



We use the WSBETA code [r2~| for calculating the eigenenergies of the WS hamiltonian 
as well as the quadrupole moments of the eigenstates. We have several parameters in the 



WS potential [O such as the depth, radius and diffuseness of the central potential and 



those of the spin-orbit potential. Among several parameter sets of the WS potential |T2 



we chose the "universal" parameter set [fL2{| . The other WS parameter sets such as the 
Blomqvist-Wahlborn, Chepurnov, Rost, and "optimal" parameter sets [0 give almost the 
same results as the "universal" parameter set, although the details of these results are not 
shown in this paper. Thus, the major results of this paper on the Ba isotopes do not change 
among those WS parameter sets. 

We assume the pairing force for the residual interaction. As our concern is the proton 
quadrupole moment, we consider only proton system hereafter. The pair potential acting 
on the individual nucleons takes up the form of 
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V pair = -A XX^yM + a{v)a{v)) (6) 

where a^f) {o,(v)) is the creation (annihilation) operator of a proton in the eigenstate v 
of the WS hamiltonian, a)(v) (a(u)) is the corresponding time-reversed operator, and the 
constant A represents the pairing gap whose value is fixed by A = G(P Jf ) in terms of coupling 
constant G and pair moment Pt defined by 

pt = £y(*) a t(i/). (?) 

u>0 



Using the Bogoliubov- Valatin transformation [|1| , we can diagonalize the hamiltonian H' = 
Pws + Vpair — XN P , where Pws is the WS hamiltonian, N p the number operator for the 
protons, and A the chemical potential. Parameters A and A are fixed by the self-consistency 
between A and P' 

A = G^u(u)v(u), (8a) 

and by the constraint on the expectation value of the particle number, Z = (N p ): 

Z = 25>(i/) 2 . (8b) 

u>0 



In Eq. (§), u and v denote the coefficients of the Bogoliubov- Valatin transformation [15 
given by 



u(u) = 2-^(l + e -^] 1/ \ (9a) 



E{u) 



e{v) - y 



with 



where eiv) stands for the single particle energy of the eigenstate v of the WS hamiltonian. 
We determined the G value for 122 Ba by the condition that the pairing gap for the 



^) = ((e(z,)-A) 2 + A 2 ) 1/2 , (10) 



protons should equal the empirical value given by JT6|: A = 12/V A ~ 1.1 MeV at /?ws=0.3, 
which /?ws value roughly corresponds to the deformation of Ba with iV = 66 as will be 
shown below. The resultant value is G = 0.16 MeV. In the summations with respect to v in 
Eq. (H), we adopted all the WS eigenstates with their eigenenergies below the barrier height 
fT2|| . We note here that the value for G depends on the number of v in these summations. In 
this paper, we assume G being independent of N, for simplicity. We have no free parameter 
except for /3ws ; i- e -> the quadrupole deformation parameter of the WS potential. The value 
for /?ws will be determined later by the experimental values of the intrinsic Qq moment. 
Refer to Fig. ||. 
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In terms of intrinsic quadrupole moment Qq for protons, we can approximately express 
B(E2:0^ -> 2t) H as 



B(E2 : 0+ - 2+) = ^OS, (11) 
with e being the proton charge. This is a good approximation for the axially symmetric 



well deformed nuclei [IB]. The quadrupole moment is given in terms of the v coefficient of 
Eq. ©by 

Q = 2 5>(z/Mz/) 2 , (12) 
where q{y) represents the quadrupole moment of the WS eigenstate v. 



q(v) = (u\^ 1 *Y 3J0 \v). (13) 

We can thus compare directly the data with the results obtained by the WS hamiltonian 
plus the BCS method. 

The intrinsic Q moment is plotted as a function of /?ws m Fig. ||. The solid curve 
represents the quadrupole moment calculated using Eq. ([12]) , which we call Q(WS+BCS). 
The filled circles denote the quadrupole moment calculated without the pairing interaction, 
which we call Q(WS). 

In general, the quadrupole moment Q(WS) increases smoothly with /3 WS because each 
occupied orbit of the ground state gains the quadrupole moment gradually with /?ws- ^ 
also happens that Q(WS) increases suddenly at a certain value of /3ws owing to a change in 



the ground state configuration. In the plot such as the Nilsson diagram [ EBJ , where the level 
energies are plotted as a function of the potential deformation, the change in the ground 
state configuration can be expected to occur at level crossings of occupied and unoccupied 
levels at the current chemical potential. 

It is clearly seen from Fig. || that Q(WS) jumps twice at /?ws ~0.025 and /?ws ~0.30. 
These jumps correspond to the level crossing where the wavefunction of the ground state 
changes its configuration. The former jump is due to the level crossing between the single 
particle orbits originating from the spherical sub-shells of M5/2 and 0^7/2. This jump is 
unimportant because the jump in Qq is too small to survive after the pairing force is switched 
on as is seen from the curve Q(WS+BCS). The jump at /3ws ~ 0.3 in Q(WS) is much bigger 
than the former one, and is due to the level crossing between the levels originating from the 
unique parity sub-shells, 0g Q / 2 and 0h u / 2 , which correspond to the levels with the Nilsson 



"asymptotic" quantum numbers (Ref. JT6J) {l[N t n z A} = | [404] and | [550], respectively, 



where is the absolute value of the z component of the total angular momentum, N t the 
total number of quanta of the deformed harmonic oscillator, n z the number of quanta in 
the oscillation along the symmetry axis of the same oscillator, and A the component of the 
orbital angular momentum along the symmetry axis [^BJ. 

Hereafter, we denote these two orbits as v\ = | + [404] and v% — | [550], respectively. 
At /9ws = 0.3, two protons, which occupy the ui and v\ orbits for /?ws < 0.3, move to the 
i>2 and z/ 2 orbits which are unoccupied by protons for /3ws < 0.3. By this movement of two 
protons, Q(WS) increases sharply at /?ws = 0.3 by 
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AQq = 2 (q(u 2 ) - q{v\)) = 1.12(b) 



(14) 



with q{vi) = —0.16 (b) and g(z^) = 0.40(b). This magnitude of AQ amounts to as much 
as 20% of the experimental value Qo=5.2 -j- 0.2 (b) for 122 Ba ||. The v\ {y 2 ) orbit has the 
largest magnitude of q in Eq. ( |i~3|) with a negative (positive) sign among the levels with 
the same sign of q near and below the Fermi surface. Therefore, emptying the v\ orbit and 
filling the v% orbit strongly affect the Qo value. 

We now consider the effect of the pairing force on Qo • The current pairing force strength 
G yields the pairing gap A = 1.34, 1.25, 0.99, 1.09, 1.06 (MeV) for /3 W s=0.0, 0.1, 0.2, 0.3, 0.4, 
respectively. Due to the level crossing between v\ and u 2 , A increases as /?ws approaches 
0.3. Owing to the pairing interaction, the sharp jumps existing in Q(WS) get smoother in 
Q(WS+BCS). The effect of the pairing force on Q depends on /3ws- The force reduces Qo 
for /3ws < 0.2 whereas it enhances Qo for 0.2 < /3ws < 0.3. The pairwise proton transfer 
takes place from Q = 9/2(g 9 / 2 ) to Q — l/2(h u / 2 ), and the probability of the pairwise transfer 
increases as /?ws approaches the level crossing point at /9ws=0.3. The quadrupole moment 
Q(WS+BCS) thus gets bigger than Q(WS) for 0.2 < (3 WS < 0.3. 

Eight horizontal lines are drawn in Fig. |2], which correspond to the experimental Qo 
moments for the Ba isotopes with N from 66 to 80. Namely, the height of each line is equal 
to the experimental value of the Qo moment. Thus, from the intersection point between each 
horizontal line and the solid curve Q(WS+BCS), we can determine /?ws corresponding to the 
experimental Qo moment for the individual isotope. It then turns out that roughly 15% of 
B(E2:0^ — > 2f) for Ba with iV = 66 is due to the pairwise proton transfer from Q = 9/2(g 9/ / 2 ) 
to Q = l/2(/zn/ 2 ). These results are not sensitive to choices of the WS parameter sets. When 



we used the WS parameter sets in [0 other than the current "universal" set, we got the 
values ranging from 15% to 20%. 

The Ba isotopes can be classified into three groups by a relation between Q(WS+BCS) 
and Q(WS): 

(a) Q(WS+BCS) < Q(WS) for iV= 76 ~ 80 

(b) Q(WS+BCS) ~ Q(WS) for iV= 70 ~ 74 

(c) Q(WS+BCS) > Q(WS) for iV= 66 ~ 68 

We should here note that A for (b) is as big as 1 MeV although Q(WS+BCS) ~ Q(WS). 

A comparison between Figs. |1] and ^| leads us to a correspondence between the results of 
IBM-2 and of the WS hamiltonian plus the BCS method: 

(a) IBM-2 with the Pauli factor gives B(E2) closer to the experimental data than the 
IBM-2 without the Pauli factor. 

(b) IBM-2 without the Pauli factor gives B(E2) closer to the experimental data than the 
IBM-2 with the Pauli factor. 

(c) Both IBM-2 with and without the Pauli factor underestimate B(E2). 

The OAI mapping used in the derivation of the IBM-2 E2 operator in Eqs. (0-^|) with 
the Pauli factor is based on the seniority scheme | JL0[| . This scheme is suitable for the case in 
which the single particle levels are nearly degenerate regarding energy at zero deformation. 
The OAI mapping thus does not take account of the intruder orbit Q = 9/2 + ((7g/ 2 ) which 
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intrudes the Z = 50 — 82 shell originating at zero deformation from below the energy gap at 
Z = 50. It is thus reasonable that the IBM-2 predictions underestimate B(E2:0 ] f — > 2~t) at 
N = 66, 68 while the IBM-2 predictions with the Pauli factor are in good agreement with 
the experimental data for the less deformed Ba isotopes with N > 76. 

We should mention here the Xe (Z = 54) and Ce [Z = 58) isotopes, though the results 
for Xe and Ce are not shown in detail in this paper. We have applied the WS hamiltonian 
plus the BCS method to the Xe and Ce isotopes in a way similar to the procedures used for 
Ba. It then turns out that, if we chose the Chepurnov set (the "universal" set) in the 
contribution of the level crossing between the orbits originating from Q = 9/2(gg/ 2 ) and the 
50-82 shell to B(E2:0^ -> 2^) for Xe with A = 66 is roughly 20% (10%). The corresponding 
quantity for Ba with A = 66 is roughly 20% (15%). For Ce with iV = 66, the effect of 
the level crossings on B(E2:0^ — > 2f) is smaller than 5% when we ignore the error in the 
experimental data. 

The B(E2:0 ] f — > 2f) value predicted by the Grodzins rule in Eq. (JI|) is only 75% of the 
experimental data for Xe with A = 66, while the prediction for Ce with A = 66 lies within 
the error bar of the experimental B(E2:0^ -> 2f ) data [§. 

The B(E2:0 ] f — > 2f) values calculated with IBM-2 without the Pauli factor are nearly the 
same as the experimental data for Xe with A = 64 — 68 M and Ce with N = 66 — 74, namely 
for the isotopes with relatively large B(E2) values. On the other hand, for the less deformed 
nuclei such as Xe with A = 70 — 80 and Ce with A = 76, IBM-2 with the Pauli factor gives 
the B(E2:0^ — > 2f) values closer to the data than IBM-2 without the Pauli factor. These 
trends are thus consistent with the IBM-2 results for the Ba isotopes in Fig. [|. This fact 
suggests that some unknown mechanism cancels the effect of the Pauli factor in Eqs. (0-0), 
and this plays a role mainly for the well deformed nuclei. To make this mechanism clear is 
beyond the scope of the current paper. We should investigate this problem further in the 
future. 

V. SUMMARY 

By the use of the Woods-Saxon hamiltonian plus the BCS method, we calculated the 
intrinsic Qo moment as a function of the deformation parameter of the WS potential. Com- 
paring the calculated Qo moment with the experimental data, we can determine the value of 
/?ws for each Ba isotope. It leads us to find a detailed way how protons occupy the deformed 
single particle orbits. By this procedure, we found that as much as 15% of B(E2:0 ] f — > 2f) 
for 122 Ba is due to the pairwise transfer from the proton intruder orbit originating from 
Q = 9/2 + (g 9 / 2 ) to the proton unique-parity orbit originating from Q = l/2~(/i n / 2 ). We 
have shown that the contribution to Qo from this pairwise transfer is roughly equal to the 
deviation of the experimental data from the theoretical predictions by IBM-2 without the 
Pauli factor. Similar deviations are found for the predictions of the Grodzins rule. We also 
found that this contribution to B(E2:0 ] f — > 2^) decreases rapidly as iV goes away from the 
mid-shell value of N = 66. This fact suggests that the pairwise transfer from the orbit with 
Q = 9/2 + (7rgg/2) to the one with Q = l/2~(7r/in/2) is responsible for the sharp increase of 
B(E2:0^ — > 2f) observed in the Ba isotope with iV = 66. 

One of the authors (K.U.) would like to acknowledge the Junior Research Associate 
Program of Japan Science and Technology Agency. 
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FIGURES 




FIG. 1. Systematics of the B(E2:0^~ — > 2+) values of the Ba isotopes. The solid squares express 
the experimental data. The open triangles present the predictions by the Grodzins rule. The solid 
(dashed) line denotes the IBM-2 results with (without) the Pauli factor. 
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Pws 

FIG. 2. Qtto -vs. /?ws plot for Ba. Proton quadrupole moment Q^q is plotted as a function of 
deformation parameter /3\ys °f the axially-symmetric deformed Woods-Saxon potential. The solid 
curve denotes <57ro(WS+BCS) which includes the pairing correlation. The filled circles express 
Qti-o(WS), the quadrupole moment calculated with the vanishing pairing gap A n = 0. Eight 
horizontal lines in the figure correspond to the experimental data. The height of each horizontal 
line is equal to the experimental quadrupole moment for each Ba isotope with the neutron number 
iV ranging from N = 66 to 80. For the sake of clarification, we have not included the errors of the 
data. The experimental values for Q w0 are 5.2 ± 0.2, 4.6 ± 0.1, 4.1 ± 0.1, 3.8 ± 0.1, 3.6 ± 0.2, 2.9 
± 0.1, 2.62 ± 0.02, 2.00 ± 0.01 (b) for 66, 68, . . ., 80, respectively [§. 
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